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1  DESCRIPTION  OF  PROGRESS 


This  report  summarizes  major  accomplishments  covering  the  period  15  March  to  13  June.  A 
paper  entitled  'A’elocity  Field  Characteristics  of  III-V  Semiconductor  Alloys"  which  was  submitted 
to  the  Journal  of  Applied  Physics  is  attached  as  Appendix  A. 

Our  aim  in  this  period  has  been  to  expand  the  work  reported  in  Appendix  .^to  study  the 
effect  of  ionized  impurity  scattering  and  strain  scattering  on  the  electron  mobility.  The  modified 
electron  mobility  will,  in  turn,  affect  velocity-field  (v-E)  characteristics  in  semiconductors  and  their 
alloys. 


Ionized  impurity  scattering  is  treated  in  a  second  order  perturbation  theory.  Our  improve¬ 
ment  over  previous  theories  is  to  convert  the  usual  integration  over  a  parabolic  energy  band  to  a 
general  Brillioun  Zone  (BZ)  integration.  This  conversion  allows  us  to  incorporate  the  effect  of  real¬ 
istic  band  structures  into  this  transport  calculation.  _ 

We  modeled  the  strain  scattering  so  the  scattering  parameters  are  calculated  from  the  empiri¬ 
cal  pseudopotential  form  factors  used  in  the  band  structure  calculation.  The  improvement  over 
previous  models  stems  from  the  universal  nature  of  the  scattering  parameters. 


The  \'CA  band  structures  have  been  used  to  obtain  group  velocity  as  a  function  of  momen¬ 
tum  and  crystal  orientation  in  various  alloys.^ The  group  velocity  curve  can  be  used  to  judge  the 
optimum  electron  energy  and  momentum  to  inject  into  ballistic  devices.  Moreover,  the  net  scatter¬ 
ing  time  due  to  all  mechanisms  can  be  calculated  as  a  function  of  electron  energy.  The  energy 
dependent  mean-free  path  is  the  product  of  the  group  velocity  and  these  scattering  times.  These 
mean-free  paths  set  an  upper  limit  on  the  length  of  the  active  region  in  a  ballistic  transport  device 
for  a  given  orientation. 

Results: 

(1)  The  effect  of  band  structure  and  the  more  accurate  BZ  integration  are  found  to  be  sub¬ 
stantial  electron  densities,  n,  in  excess  of  10^^  cm"^.  At  n  =  10^^  cm'^,  band  struc¬ 
ture  effects  alone  reduce  the  electron  mobility  by  about  25%  relative  to  those  calculated 
in  the  constant  effective  mass  approximation  for  an  electron  temperature  of  1000  K. 
This  effect  needs  to  be  included  in  any  accurate  calculation  of  v-E  behavior. 

(2)  Calculations  of  group  velocity  from  VGA  band  structure  suggest  InAs  or  InSb  com¬ 
pounds  will  be  better  suited  to  ballistic  electron  devices  than  GaAs.  As  seen  from  the 
enclosed  curves,  the  maximum  group  velocity  always  occurs  in  the  (100)  directions,  and 
ImAs  is  about  70%  larger  than  that  in  Ga-As.  Recent  experiments  [Levi,  et  al.,  Appl. 


I’hys.  Lett.,  (1986)  to  be  publislied;  also  find  that  liiAs-based  devices  are  .about  7  times 
more  elficient  than  those  based  on  Ga^-Vs.  We  still  must  calculate  these  elTocts  including 
the  full  CP.A  band  structures.  Our  calculations  further  show  that  in  (100)  is  larger 
than  that  in  cither  ( 1 10)  or  (111). 

Because  lime  is  now  available,  we  are  now  engaged  in  calculation  of  CPA  band  struc¬ 

tures  of  various  alloys.  The  effect  of  real  part  of  self  energy  will  change  the  density  of  states  and 
group  velocities  at  liiglier  energies.  These  changes  will  be  included  in  the  calculation  of  alloy, 
impurity,  strain  and  phonon  scattering.  .Vlso  these  calculations  will  include  anisotropic  scattering. 
Then,  the  peak  velocity  and  the  saturation  velocity  will  be  obtained  for  various  alloys. 

II  EQUIPME.NT  PURCIL\SED  OR  CONSTRUCTED 

.None 

III  TRIPS.  MEETINGS,  P.VPERS,  AND  VISITS 

Dr.  Arden  Slier  and  Srinivasan  Ivrishnamurthy  attended  an  .\merican  Physics  Society 
Conference  in  Las  Vegas.  Nevada,  March  31  to  April  4.  1986.  The  following  paper  has  been  sub¬ 
mitted  to  Journal  of  Applied  Physics:  'Velocity-Field  Characteristics  of  III-V  Semiconductor 
Alloys,"  by  S.  Krisiinamurthy,  A.  Sher,  and  A.-B.  Chen. 

IV  PROBLEMS  OR  AREAS  OF  CONCERN 

•None 

V  DEVIATION  FROM  PLANNED  EFFORT 

None 

\1  FISCAL  STATUS 
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The  total  contract  funding  for  the  three-year  period  is  $611,296.  Of  this,  $195,287  was  allo¬ 
cated  to  the  first  year.  Including  the  burden,  approximately  $22,000  was  intended  to  pay  the  con¬ 
sulting  fee  of  .\n-Ban  Chen  leaving  approximately  $173,000.  In  the  first  year  ending  31  May  1986, 
we  spent  approximately  $192,000.  Year  two  has  been  allocated  $203,655,  with  $20,284  intended  to 
pay  the  consulting  fee  of  .Aji-Ban  Chen. 
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ABSTR.ACT 

We  have  calculated  the  velocity-field  characteristics  of  semiconductor 
alloys  based  on  realistic  band  structures;  we  obtained  the  band  structures  and 
alloy-scattering  rates  from  a  generalization  of  the  coherent  potential  approxi¬ 
mation  method.  Although  we  use  proper  band  structures,  w’e  still  consider  a 
single  electron-temperature  model.  The  results  agree  surprisingly  well  with 
experiments,  and  suggest  that  InP-based  alloys  are  good  candidates  for  high¬ 
speed  devices. 
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I.  INTRODUCTION 


In  many  semiconductors,  the  electron  drift  velocity,  v.  increases  nonlineariy  under 
the  influence  of  high  fields  to  reach  a  peak  value,  Vp,  at  a  threshold  field.  Ej.  and  then 
decreases  with  applied  E,  giving  rise  to  a  negative  differential  resistance.  For  device 
applications,  one  must  know  the  effect  of  material  (band-structure)  parameters  on  hot 
electron  transport.  In  this  paper,  we  examine  the  influence  of  the  band-structure  effects 
on  velocity-field  characteristics  of  III-V  ternary  alloys.  Such  information,  calculated 
differently,  is  already  available  for  several  III-V  and  II-\'I  compounds  and  some 
alloys. Our  improvement  over  previous  efforts  is  to  use  accurate  band  structures  and 
alloy  scattering  rates  for  many  systems.  Our  aim  is  to  examine  systematically  the  factors 
which  limit  the  values  of  Ej  and  Vp  and  to  suggest  likely  materials  for  high-speed  device 
applications. 

n.  MODEL  FOR  v-E  CHARACTERISTICS 

For  this  comparative  study,  we  have  adopted  a  simple  hot  electron  theory"  that 
assumes  a  '“drifted  Maxwellian  ’  electron  distribution  with  a  single  electron  temperature. 
T,.  This  distribution  would  be  established  only  if  there  were  an  interaction,  like 
electron-electron  scattering,  that  drove  the  electron  system  into  internal  equilibrium  at  a 
rate  fast  compared  to  the  rate  of  energy  exchange  with  the  heat  bath.'*  No  such  inter¬ 
action  has  ever  been  identified  because  it  requires  very  high  electron  density.  Neverthe¬ 
less.  the  predictions  of  the  drifted  Maxwellian  distribution  are  surprisingly  accurate,  and 
it  serves  as  a  good  vehicle  for  this  systematic  study,  so  we  are  adopting  it  in  this  paper. 

Under  high  field,  the  electrons  gain  energy  when  they  are  accelerated  from  the  con¬ 
duction  band  minimum  to  higher  energy  levels.  The  average  electronic  energy  in  a 
steady  state  is  often  considerably  higher  than  that  of  the  lattice  excitation  energy  per 
particle.  The  electron  temperature.  T..  which  represents  the  average  electron  energy,  is 
therefore  much  higher  than  the  lattice  temperature,  T.  For  a  direct-gap  semicontluctor. 


the  mobility  in  the  lower  valley,  /ij,  is  much  larger  than  that  in  the  higher  valleys,  /i.>,  so 
the  drift  velocity  can  be  approximated  by" 


V  =  /iiE/(l+Re 


(1) 


where  R  is  the  ratio  between  the  densities  of  states  in  the  upper  and  lower  valleys, 
respectively,  and  AE  is  the  energy  difference  between  the  minima  of  these  two  vallevs.  In 
steady  state,  the  energy  gained  by  an  electron  from  the  field  must  be  balanced  by  the 
energy  lost  to  the  lattice 


evE  =  (3/2)  k(T,  -  TI/te 


where  Te  is  the  electron-lattice  energy  relaxation  time.  Equations  (1)  and  (2)  combine  to 
yield  T.,  which  in  turn  yields  the  v-E  characteristic  curve  and  hence  the  threshold  values 
Ex  and  Vp.  These  values  have  the  following  approximate  forms:  Vp  oc  v/tiAE^^  and 
Ex  VR-^/ Mi'E' 

Alloy  effects  come  not  only  from  the  concentration  dependence  of  the  band  parame¬ 
ters  R  and  AE,  but  also  from  the  effect  of  alloy  scattering  on  the  mobility.  We  take  /i; 
to  be 


Mr‘  =  (T)  +  (T,) 


where  Hq  is  the  mobility  limited  by  al'  scattering  mechanisms  except  alloy  scattering  and 
is  taken  as  the  concentration-weighted  average  of  the  experimental  values  of  the  two  con¬ 
stituent  compounds.  is  the  alloy-scattering-limited  mobility  and  is  calculated  using  a 
generalized  Brooks'  formula.'^ 
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Finally,  we  assume  that  electrons  lose  energy  to  the  lattice  mainly  through  polar 
optical  phonon  scattering  and  calculate  from  Eq.  (3. 6. 26)  of  Ref  1.  The  electronic  and 
lattice  parameters  needed  for  this  estimate  are  all  taken  to  be  the  concentration-weighted 
averages  of  the  pure  crystal  constituents. 

ra.  BAND  STRUCTURES 

To  calculate  the  alloys’  band  structures,  we  must  begin  with  those  of  the  pure  com¬ 
pounds.  We  have  demonstrated^  that  high-quality  band  structures  for  [I-\  I  and  III-\ 
compounds  can  be  obtained  using  a  minimum  set  of  four  sp^  Slater  orbitals  per  atom  in 
a  semi-empirical  calculation.  First,  the  empirical  pseudo-potential  form  factors  are  used 
to  construct  a  tight-binding  (TB)  Hamiltonian  II  in  this  minimum  set.  H  is  then 
transformed  into  a  zeroth  order  Hy  in  an  orthonormal  basis.  Then  a  perturbative  Hamil¬ 
tonian  having  a  first-neighbor  TB  form  and  a  site-diagonal  spin-orbit  Hamiltonian  are 
added  to  obtain  a  fine-tuned  band  structure.  This  band  structure  accurately  reproduces 
the  energy  separation  between  the  lowest  and  the  next  lowest  minima  in  the  conduction 
bands,  the  band  gap  Eg,  and  the  conduction  band  effective  mass  m*. 

Figure  1  shows  an  example  of  a  band  structure  for  InAs.  The  quality  of  the  band 
structure  can  be  seen  to  be  comparable  to  the  best  ones  available.'  The  band  structure 
obtained  in  this  way  is  superior  to  the  results  of  the  usual  first-  or  second-neighbor  TB 
approach  because  our  method  includes  all  the  long-range  intera!  tions.  Figures  2(a)  and 
2(b)  show  the  lowest  conduction  band  in  FX  and  FL  directions  for  In.\s.  InSb.  InP  and 
Ga.\s.  GaSb  compounds,  respectively.  The  phenomena  we  are  examining  are  only  sensi¬ 
tive  to  energy  variation  of  states  relative  to  the  conduction  minima,  so  all  the  curves  are 
drawn  with  their  conduction  band  edges  aligned.  Because  the  InP,  InSb  conduction 
bands  cross  in  Figure  2(a),  their  alloys  will  have  larger  /t,  but  smaller  AE  than  does  InP. 
The  peak  velocity  in  these  alloys  can  be  larger  than  that  of  InP.  However,  both  /tj  and 
AE  will  be  larger  in  InP.Ajs  alloys  than  that  in  InP.  One  can  expect  a  moderate  increase 
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in  and  Ey  '  Because  //j  and  AE  are  almosc  the  same  in  GaAs  and  GaSb.  electron- 
phonon  interactions  then  determine  the  change  in  v-E  characteristics. 

Table  I  lists  the  important  band  parameters  of  the  six  pure  compounds 

included  in  this  study.  Also  included  is  the  minimum  energy  a  hot  electron  needs  to 
create  an  electron-hole  pair.*  Eg  (I+'2Q)  (l+a)',  where  a  is  the  ratio  of  the  electron  mass 
to  the  heavy-hole  mass. 

These  pure-crystal  band  structures  are  the  input  to  the  alloy  electronic  structure 
calcu'ation  using  the  coherent  potential  approximation  (CP.\).^  Both  diagonal  disorder 
due  to  atomic  term  value  differences  and  off-diagonal  disorder  due  to  bond-length 
fluctuations  are  treated.  The  method  is  a  generalization  of  the  tight-binding  ''molecular 
CP.\  "  of  Hass  et  al.®  to  our  better  band  structures.  From  the  CP.\  self-energy,  one  can 
obtain  the  alloy  band  gap.  m*,  AE,  R.  and  /aj  as  a  function  o'"  the  alloy  concentration  x. 
These  quantities  are  then  used  to  calculate  the  v-E  characteristics  as  we  described  in 
Section  11. 

T\\  EX.AMPLE:  GaJni_^As 

As  an  example,  we  present  the  results  of  our  calculations  for  Ga^In]_x.\s  in  some 
detail  in  this  section.  Figure  3  plots  AE,  Ep,  and  Vp  at  two  temperatures  as  a  function 
of  X.  If  AE  is  larger  than  Ep,  avalanche  breakdown  is  expected  to  occur  before  Vp  is 
reached.  Thus,  the  velocity-field  plot  near  Vp  is  only  meaningful  for  the  range  of  x  where 
AE  <T  Ep.  The  maximum  value  of  Vp  that  can  be  obtained  for  this  alloy  is  at  x  —  0.-47. 
The  peak  vplocity.  Vp,  is  seen  to  be  a  nonlinear  function  of  x.  Experimental**  values  of 
Vp  (open  circle)  are  also  plotted. 

Because  we  set  the  mobility  of  the  satellite  valley  equal  to  zero  in  Eq.  (1).  we 
should  not  expect  the  model  to  be  accurate  for  electric  fields  much  higher  than  Ey 
Nevertheless,  it  is  interesting  to  examine  the  v-E  curves  produced  by  the  model  for 
different  x  values  shown  in  Fig  4.  Because  the  mobility  of  InAs  is  larger  than  that  of 
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GaAs.  the  slope  in  the  low  field  region  decreases  with  x.  For  the  same  reason.  Ej 
increases  with  x.  This  prediction  of  a  larger  Vp  value  and  a  smaller  Ej  at  lower  GaAs 
content  has  been  observed  experimentally.**  At  x  =  0A7.  the  Vp  value  is  about  o0% 
larger  than  that  of  GaAs,  once  again  in  agreement  with  experiments.'**'*' 

The  effect  of  the  lattice  temperature  on  the  v-E  characteristics  is  illustrated  in 
Fig.  .5.  The  curve  sharpens  at  low  T  and  the  Vp  value  increases,  because  both  //g  and 
are  larger  at  lower  T. 

Figure  6  depicts  two  effects:  the  influence  of  alloy  scattering  on  the  v-E  charac¬ 
teristics  and  the  errors  caused  by  the  effective  mass  approximation  (ENL\).  All  three 
curves  are  for  x  =  0.5  and  at  T  =  300  K.  Curve  (a)  does  not  include  alloy  scattering 
and  curve  (c)  does,  but  both  are  based  on  EMA.  The  presence  of  alloy  scattering  reduces 
the  effective  mobility,  so  curve  (c)  starts  with  a  smaller  slooc  and  has  a  smaller  Vp  and  a 
long  tail  as  it  approaches  Vp,  causing  its  Ej  to  be  larger  than  curve  (a).  The  reduction  of 
Vp  by  alloy  scattering  in  this  case  is  about  30Gc.  Curve  (b)  represents  the  results  of  a 
more  detailed  calculation  using  the  actual  density  of  states  (DOS).  Instead  of  Eq.  (1),  the 
drift  velocity  in  this  calculation  is  obtained  from  v  =  i7E,'(ne),  where  n  is  the  electron 
density  and  a  is  given  by 


(7  =  (•2,3)e-  /  tr{€)T(i)  (-d('d€)de  ,  (3) 

with  r(e)  being  the  average  momentum  collision  time,  f  the  electron  distribution,  and 
the  energy-shell  average  of  the  square  of  the  group  velocity 

=  E  lE.{k)  1=  (-1) 

nk 
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Thus  curve  (b)  does  not  assume  a  zero  mobility  for  the  liigher  valleys  as  was  assumed  in 
curves  i^a)  and  fc),  W’e  therefore  infer  that  the  effective  mobility  in  curve  (b)  lies  between 
those  of  curves  (a)  and  (c).  Comparison  between  curves  (b)  and  (c)  shows  that  the  more 
elaborate  calculation  tends  to  have  a  moderate  increase  in  Vp  and  a  small  reduction  in 
E-j".  largely  because  the  bands  increase  with  k  more  rapidly  than  quadraticallv  before 
they  bend  over  toward  the  satellite  valley.  Thus  as  a  function  of  k.  the  effective  mass 
actually  first  becomes  smaller,  then  larger  until  the  inflection  point  is  reached,  where  the 
curvature  changes  sign  and  the  effective  mass  becomes  negative.  The  alloy  mobility 
varies  roughly  as  (m*)”'^'",  so  the  net  effect  of  the  smaller  m  in  the  critical  region  is  to 
produce  curve  (b).  Moreover,  the  real  part  of  the  alloy  scattering  self-energy  must  also 
play  a  role,  because  it  too  modifies  the  energy  bands.  Thus.  Monte  Carlo  calculations, 
which  treat  scattering  on  fixed  energy  band  structures,  cannot  be  expected  to  reproduce 
the  same  numerical  result. 

V.  OTHER  TERNARY  ALLOYS 

We  have  carried  out  the  above  calculations  for  the  direct-gap  ternary  alloys 
GalnSb,  In.\sP.  InAsSb,  GaAsSb,  and  InPSb.  The  results  are  summarized  in  Fig.  7.  the 
quantities  plotted  are  the  values  of  Vp  and  Ex~^  which  are  respectively  the  Vp  and  Ej”' 
normalized  to  GaAs  values.  Also  plotted  are  the  values  of  AE  and  Ep.  Because  only  a 
few  alloys  satisfy  the  condition  Ep  >  AE.  the  peak  values,  Vp,  shown  can  actually  be 
reached  before  avalanche  breakdown  takes  place.  The  circles  are  the  available  experi¬ 
mental  values"'^^''^  for  Vp.  Considering  the  simplifications  involved  in  the  model,  the 
agreement  with  experiment  is  remarkable. 

It  is  interesting  to  note  that  Vp  is  greater  in  InPSb  than  in  pure  InP  and  InSb. 
although  alloy  scattering  tends  to  reduce  Vp.  Similar  results  are  also  seen  in  InP.\s  alloy. 
These  results  can  be  explained  as  follows:  the  phonon  scattering  rate  in  InP  is  larger 
than  that  in  either  InSb  or  InAs  However,  because  of  its  larger  effective  mass,  the 
mobility  in  InP  is  much  smaller  than  those  in  InSb  or  InAs.  The  combined  effect  gives 


rise  to  a  larger  Vp  for  the  alloys.  The  ratio  of  the  peak  velocity  to  the  energy  dissipated 
by  an  electron  traveling  at  this  velocity  is  Ej'*,  and  so  E-j-”'  measures  the  efficiency  of 
an  alloy  relative  to  GaAs  in  high-speed  device  applications.  Our  calculation  of  E-j-”*  also 
agrees  very  well  with  the  available  experimental  data. 

VI.  SUMMARY 

We  have  carried  out  a  comparative  evaluation  of  the  influence  of  band  structure 
fe^'  -“s  of  several  III-V  compound  alloys  on  velocity-field  characteristics  using  detailed 
band  structures  and  alloy  scattering  rates.  Our  results  suggest  that  InP-based  alloys  are 
even  more  promising  candidates  for  high-speed  devices  than  Ga.\s.  A  more  rigorous 
study  is  in  progress  to  improve  the  accuracy  of  our  calculations  and  further  check  this 
prediction. 
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FIGURE  CAPTIONS 


1.  Band  structure  of  InAs.  Energy  in  e^  . 

2.  Lowest  conduction  band  of  (a)  InP.  In.\s,  InSb  and  (b)  GaAs.  GaSb  compounds.  All 

energies  measured  in  eV  with  respect  to  their  respective  conduction  band  minima. 

3.  \’ariation  of  AE,  E^,  v.  as  a  function  of  Ga.\s  concentration.  \  in  Ga-^In|_^A.s  alloys. 

Experimental  values  (open  circle)  of  Vp  are  from  Ref  11. 

4.  \'elocity-field  characteristics  of  Ga^Ini_xAs  alloys. 

■5.  Temperature  dependence  of  v-E  characteristics  of  Gao  ^Iiif)  ;^.\s 

6.  Eflect  of  band  structure  and  alloy  scattering  on  v-E  characteristics.  Results  ol  the 

calculations  (a)  in  ENL-V  without  alloy  scattering,  (b)  in  EM.\  with  alloy  scattering, 
and  (c)  with  exact  DOS  and  alloy  scattering. 

7.  \’ariations  of  AE  (solid  line — in  lower  half),  Ep  (dotted  line — in  lower  half).  Vp  (solid 

line— in  upper  half),  and  Ex”'  (dotted  line — in  upper  half)  are  shown,  Vp,  Et”'  are 
peak  velocity  and  threshold  field  in  the  unit  of  respective  GaAs  values.  The  experi¬ 
mental  values  of  peak  velocity  (open  circle)  and  threshold  fields  (filled  circle)  are 
also  shown. 
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Table  I,  Important  band  parameters  of  compounds. 


Parameter 

GaAs 

In  As 

InP 

GaSb 

InSb 

EgleV) 

1.52 

0.42 

1.42 

0.S6 

0.25 

m*(  m^) 

0.067 

0.026 

0.072 

0.0.50 

0.018 

AE(eV) 

0.2S2 

l.Ii 

0.S2 

0.27 

0.73 

R 

42.95 

186.02 

45.02 

200.10 

1493.14 

Ep(eVI 

1.70 

0.44 

1.54 

0.98 

1.26 

